Local Spin Susceptibility of the 5=1/2 Kagome Lattice in ZnCu3(OD)6Cl2 
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We report single-crystal NMR investigation ol the nearly ideal spin S = 1/2 kagome lattice 
ZnCu3(OD)6Cl2. We successfully identify NMR signals originating from the nearest-neighbors 
of Cu^+ defects occupying Zn sites. From the Knight shift measurements, we demonstrate that 
weakly interacting Cu^'*' spins at these defects cause the large Curie- Weiss enhancement toward T = 
commonly observed in the bulk susceptibility data. We estimate the intrinsic spin susceptibility 
of the kagome planes by subtracting defect contributions, and explore several scenarios. 



Persistent search for frustrated spin-liquid systems 
without a magnetically ordered ground state has re- 
cently resulted in the discovery of several model materials 
with spin S = 1/2, including the kagome [3H5], triangular 
[5], and FCC [3 [S] lattice systems. Intense experimen- 
tal and theoretical efforts are underway to understand 
their ground state properties, and to identify potential 
signatures of exotic spin-liquid states. Among the most 
promising candidates of spin-liquids is herbertsmithite 
ZnCu3(OH)6Cl2, which contains planes of a kagome lat- 
tice formed by Cu^^ ions with S = 1/2 (see Fig. 1 for 
the crystal structure) [31 SI [5HT^ . The antiferromagnetic 
Cu-Cu exchange interaction J is estimated to be as large 
as J/ks ~ 190 K [4 , yet there is a firm consensus that 
no hint of magnetic long range order exists at least down 
to 50 mK Hin]. 

Nonetheless, the mechanism behind the realization of 
such a non-magnetic ground state in ZnCu3(OH)gCl2 has 
been controversial (see [TS] for a recent review). The 
root cause of the controversy stems from the difficulty 
in pinpointing the location of defects, and understand- 
ing their potential influence on various physical proper- 
ties. For example, magnetic Cu^+ defects occupying the 
non-magnetic Zn^"*" sites outside the kagome plane might 
induce a large Curie- Weiss contribution xcw to the bulk- 
averaged spin susceptibility Xbuik, and mask its intrinsic 
behavior. In fact, Xbuik always exhibits a large enhance- 
ment below about 50 K, as shown in Fig. 2a. The ob- 
served behavior is different from the predicted behavior 
of Xkagome for ideal kagome antiferromagnets with van- 
ishing spin susceptibility at T = [TBl - B^ . Furthermore, 
if non-magnetic Zn^+ defects occupy the magnetic Cu^"*" 
sites within the kagome plane, spin singlets may be in- 
duced in their vicinity toward T — Q jT^] , and help realize 
the non-magnetic ground state. These uncertainties over 
the nature of defects prevented extraction of the intrinsic 
Xkagome from Xbulk ■ It is therefore necessary to untangle 
the effects of defects based on real space probes. 

In this Rapid Communication, we take advantage of 
recent successful growth of high quality single crystals 
of deuterized ZnCu3(OD)6Cl2 [H], and report a single- 
crystal NMR investigation. Thanks to very sharp NMR 
lines in a 45 mg single crystal, we identified NMR sig- 
nals originating from the 6 nearest-neighbor (n.n.) and 



12 next-nearest-neighbor (n.n.n.) sites of Cu^+ de- 
fects occupying the Zn sites. From their relative inten- 
sities, we show that magnetic Cu^+ defects replace non- 
magnetic Zn^"*" with a probability of 14 ± 2%. A Curie- 
Weiss behavior of the NMR Knight shift observed at 
the n.n. sites indicates that the large Curie-Weiss en- 
hancement commonly observed for Xbuik below ^ 50 K 
arises from these defects. By subtracting the contribu- 
tion of the Curie- Weiss contribution xcw from Xbuik^ we 
explore the intrinsic behavior of spin susceptibility of the 
kagome planes, Xkagome- 

In Fig. 3, we present representative NMR line- 
shapes observed in an external magnetic field of B^xt = 
8.4 Tesla applied along the c-axis, i.e. normal to the 
kagome plane. The nucleus has a nuclear spin 1 = 1, 
and hence we would expect to observe a pair of NMR 
transitions between the nuclear spin state = ±1 
and Iz = 0, if there was no defect. In general, one 
can express the resonant frequency of each site as 
f± = jnBtotai ± Vq, where the nuclear gyromagnetic 
ratio is 7n/27r = 6.53566 MHz/Tesla, and Btotai repre- 
sents the total magnetic field seen by each nuclear spin. 
The second term, ±i/q, represents contributions of the c- 
axis component of the nuclear quadrupole interaction for 
the aforementioned = ±1 to = transitions. The 




FIG. 1. (Color Online) Left: local configuration of ions sur- 
rounding the Zn^^ site in ZnCu3(OD)6Cl2. Right: the c-axis 
view of the kagome plane. 
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Uq is proportional to the second derivative of the total 
Coulomb potential at the sites, and small differences 
in the local structure lead to distinct values of the NMR 
line splitting Sj/q. 

In reality, we observe three pairs of NMR signals, 
A±, B±, and C±. This means that there are three dis- 
tinct sites in ZnCu3(OD)6Cl2 with locally different 
environment. Based on the relative intensity of well- 
resolved peaks A+, B_|_, and C+ at 295 K, we estimate 
the population of A, B, and C sites as 14 ± 2%, 28 ± 4%, 
and 58 ±4%, respectively. From the line splitting, we ob- 
tain = (/+ - /_)/2 = 48, 62, and 45 kHz for A±, B±, 
and C± peaks, respectively. Within experimental uncer- 
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FIG. 2. (Color Online) (a) (+): The bulk averaged mag- 
netic susceptibility xtuik measured by SQUID in 0.2 Tesia for 
a collection of small single crystals with random orientations. 
The solid curve represents extrapolation of the Curie- Weiss 
contribution xcw with 9buik = —1.2 K estimated by fitting 
Xbuik below 3 K. (•): NMR Knight shift measured 
at the A site, which exhibits nearly identical temperature 
dependence as Xcw- (°): The intrinsic spin susceptibility 
within the kagome plane, Xkagome = Xbuik - Xcw- (b) (o): 
the same Xkagome plotted in (a). The solid curve through the 
data points represents the best fit with a small gap (see main 
text). The dotted curve shows the theoretical prediction [18) 
with J/ks = 180 K. (•): Xkagome estimated from the best 
fit of Xb-uik below 8 K to a summation of the T-linear term 
from Dirac Fermions and a Curie- Weiss term. (A): an alter- 
native scenario of Xkagome with a finite gapless excitations c 
at T = 0. 



tainties, we didn't observe any noticeable temperature 
dependence in Vq. Thus the temperature dependence of 
f± arises from magnetic effects through ^nBtotai- 

The B total at each site is a summation of the ex- 
ternal field Bext — 8.4 Tesla and the hyperfine magnetic 
field Bhf from nearby electron spins polarized by B^xt^ 
i.e. Btotai — Bext + Bfif. Since B^f is proportional to lo- 
cal spin susceptibility Xspin in the vicinity of the observed 
site, one can gain insight into the position dependent 
Xspin within the kagome plane by keeping track with the 
temperature dependence of f± at different sites. A strik- 
ing feature of Fig. 3 is that A± peaks progressively shift 
to lower frequencies, and split off below 10 K. As shown 
below, this is because A sites are under the direct influ- 
ence of the Cu^+ defect spins at n.n. Zn sites. 

In Fig. 4, we summarize the temperature dependence 
of the 2d NMR Knight shift = BufJ B,xt = {Btotai - 
Bext)/ Bext, which is related to local spin susceptibility as 



hfxipin+Kchem- A)^'f IS the hyperfine coupling 



constant between the observed nuclear spin and the 
electron spin at the ith Cu site in their vicinity, and xipm 

is the local spin susceptibility at the ith Cu site. A^j!^j is 
negative in the present case, hence we reverse the vertical 
axis of Fig. 4. Kohem is a small temperature independent 
chemical shift originating from the orbital motion of elec- 
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FIG. 3. (Color OnHne) Representative NMR lineshapes 
observed in Bext = 8.4 Tesla applied along the c-axis. We 
observe three pairs of quandrupole-split NMR lines, A±, B±, 
and C± . Downward arrows mark the resonant frequencies f± 
for A±. Notice that the horizontal scale is changed for (b). 
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trons. The results in Fig. 4 indicate that the Knight shift 
'^K{A) at the A sites is anomalously large compared with 
'^K{B) and ^K{C) at B and C sites. Furthermore, the 
temperature dependence of ^K{A) obeys a single Curie- 
Weiss law in the entire temperature range with a Weiss 
temperature, 9a = —3.5 K. 

In Fig. 2, we compare the temperature dependences of 
'^K{A) and Xtuik- It is interesting to note that '^K{A) 
asymptotes to Xbuik at the low temperature limit. In 
fact, the low temperature behavior of Xbuik can be ap- 
proximated by a Curie- Weiss law xcw ~ 0.12/(T + 
(^buik) emu/mol-f.u. with a vanishingly small Weiss tem- 
perature, Obuik — —1.2 K; the extrapolation of xcw to 
higher temperature reproduces the temperature depen- 
dence of ^if(A), as shown by a solid curve. The slight 
overestimation of Oa = —3.5 K in B^xt = 8.4 Tesla by 
NMR over 8\juik = —1.2 K should be attributed to a sat- 
uration of defect spin polarization in high fields below 
T ~ HsBext/kB ^ 5.6 K. We confirmed that ^K{A) at 
4.2 K measured in Bext — 5.3 and 3.2 Tesla is indeed 
somewhat larger, and the field dependence is consistent 
with that of bulk magnetization. 

Thus we have identified the source of the large en- 
hancement of Xbuik observed below 50 K as arising from 
defect moments neighboring A sites whose behavior 
can be modeled with a Curie- Weiss law. The vanishingly 
small Weiss temperature implies that interactions among 
these defect moments are very weak. This weak coupling 
may explain the effectiveness of a simple Curie- Weiss law, 
even though a distribution of coupling strengths would 
be expected for a random arrangement of defects. The 
next question to address then, is: where are these defects 
located? We recall that A sites account for 14±2% of the 
overall NMR signals, and the intensity of B sites is 
precisely twice greater, 28 ± 4%. In view of the fact that 
each Cu^+ defect occupying the Zn sites has 6 n.n. and 
12 n.n.n. sites (see Fig. 1), we assign A and B sites as 
the n.n. and n.n.n. sites of the Cu^"*" defects, respectively. 
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FIG. 4. NMR Knight shift '^K at three distinct sites. 
Solid curve represents the best Curie- Weiss fit. 



These site assignments also imply that 14 ± 2% of the Zn 
sites in our single crystal are occupied by weakly inter- 
acting Cu^"*" spins, and there is no NMR signature of 
Zn^+ anti-site defects occupying the Cu^"*" kagome sites. 
Our findings are consistent with the recent site-specific x- 
ray anomalous dispersion measurements combined with 
Rietvelt refinement, from which the actual composition of 
our sample was estimated as (Zuq ssCuq i5)Cu3(OD)6Cl2 

m- 

We gain additional insight into the nature of defects 
from the sign of the hyperfine interaction, A^f- Our 
observation, '^K{A) < 0, implies that the hyperfine 
field Bf^f at A sites, induced by the defect mo- 
ments polarized along the c-axis in B^xt, is negative. 
From the comparison of the temperature dependence of 
'^K{A) and xcw, we estimate A^f oc K) / d{xcw) = 
—24 kOe/nB- On the other hand, in the present ge- 
ometry employed for our NMR measurements, it is easy 
to show that the dipole hyperfine field from Zn sites, 
^thpoie _ _j_Q g^ kOe//iB, has a positive component along 
the c-axis. Fermi's contact hyperfine interaction induced 
by direct overlap of the Is orbital and Cu 3d^ or- 
bitals at Zn sites is also positive. These considerations 
suggest that a sizable spin polarization obeying the same 
Curie- Weiss behavior also exists in the Cu3(OD)3 trian- 
gles within the kagome planes above and below the Zn 
defect sites; the negative hyperfine field at A sites 
originates from these spin polarizations. This conclusion 
also explains the fact that 14 ± 2% of isolated Cu^+ de- 
fects alone would underestimate the magnitude of xcw 
by a factor of ~ 2. 

In contrast with the Curie- Weiss behavior of A sites, 
the Knight shift ^i^(C) at C sites is very small, and 
doesn't show a large enhancement toward T = 0, as 
shown in Fig 4. In other words, at least 58% of Cu sites 
within the kagome planes show no hint of enhancement of 
local spin susceptibility Xspin toward T = 0. This means 
that the defect moments do not significantly couple to 
the Cu ions in the kagome plane beyond the n.n.n. Cu 
sites, with which the 2D C sites are bonded. Since we 
can't resolve the B sites from C sites below~ 50 K, 
28% of Cu sites adjacent to B sites don't have a large 
enhancement of Xspim either, although we can't entirely 
rule out a mild enhancement. 

We have established that the Curie- Weiss contribution 
Xcw to Xbuik originates from the immediate vicinity of 
the Cu^"*" defect moments occupying the Zn sites. This 
finding paves a new avenue to analyze the intrinsic spin 
susceptibility within the kagome plane, Xkagome, by sub- 
tracting Xcw from Xbuik, i-e. Xkagome = Xbuik - XCW ■ 

In earlier studies of various quantum antiferromagnets, 
Xcw was generally subtracted from Xbuik to reveal in- 
trinsic behavior, e.g. [23]. Because of the uncertainties 
in the origin of defects, such analysis of Xkagome has not 
been fully explored in the present case. In Fig. 2(a), we 
present Xkagorae = Xbuik - XCW, together with Xbuik and 
Xcw- We also present the same Xkagome in Fig. 2(b) in 
an expanded scale. 
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Once we subtract xcwj a qualitatively different behav- 
ior of Xfcagome emerges; Xkagome Saturates below 100 K, 
then decreases with temperature toward T = 0. The 
temperature dependence of Xkagome is consistent with 
the theoretical prediction based on a high temperature 
series expansion with J/fcs ^ 180 K as shown by a dot- 
ted curve. Our finding near T = is also consistent 
with earlier reports based on -^^Cl and ^^O powder NMR 
that Xspin at some fraction of Cu sites decreases toward 
T = [TTJ [12] . However, dramatic NMR line broadening 
hampered the efforts to probe the bulk-averaged behavior 
below 50 K in these earlier powder NMR studies. 

The solid curve in Fig. 2b represents the best fit to 
an empirical formula Xkagome = A ■ exp{—A/kBT)/{T + 

Okagome) with a Small gap A/fcs 7 K and 9kagome ^ 

690 K. The fit is reasonably good, but we should be cau- 
tious in interpreting its implications. Since we estimated 
Xcw from Xbuik below 3 K, and subtracted the former 
from the latter, in essence, we are forcing Xkagome to 
reach zero at 3 K. In other words, the estimated values 
of Xkagome sct the lowcr bound of the intrinsic spin con- 
tributions at kagome sites. Thus we can't rule out the 
possibility that this small gap is an artifact of our anal- 
ysis. In this sense, we should consider A ^ 0.04J as the 
upper hound of any possible spin gap. 

In fact, we noticed that the temperature dependence 
of Xkagome near T = is also consistent with a T-linear 
behavior predicted for gapless excitations of the Dirac 
Fermion in algebraic spin liquids, Xkagome = AjjpT, 
where Adf = Q.&N Ai^i^B / J"^ emu/Kmol-f.u.pUl [^. It 
is worth noting that the same Dirac Fermion model suc- 
cessfully accounts for the Cspin ^ behavior of the spe- 
cific heat data below 30 K with J/fcs - 190 K [2011^1] . 
Inspired by these promising signatures, we also fitted 
Xbuik as a summation of the Dirac Fermion contributions 
Xklgome ^ ^DpT and defect contributions Xcw below 
8 K, and estimate Xkfgome = Xbuik - Xcw- turned out 
that Xcw nearly identical with the previous estimate, 
hence Xkagome shows very similar behavior, as shown in 



Fig.2b. The best fitted value, Adf = 2.0 x lO"'' emu/K 
mol-f.u., implies J/kB ~ 220 K. This is in reasonable 
agreement with J/fc^ = 180 ^ 190 K mentioned above. 

While the predictions based on Dirac Fermions with 
zero density of states at the Fermi energy account for 
both Xkagome and Cspin fairly well, we caution that the 
nature of low lying states of spin liquids is still contro- 
versial. For example, spin liquids with a spinon Fermi 
surface on triangular lattices are known to exhibit a fi- 
nite density of states [6l |25] . An analogous scenario, how- 
ever, seems highly unlikely in the present case, because 
Cspin exhibits positive curvatures in a broad tempera- 
ture range below 30 K f?' , instead of a T-linear behavior 
expected for a finite density of states. We estimate the 
upper bound of the T-linear term as Cspin ^ O.IT J/mol- 
f.u. below 1 K in 14 Tesla [1]. This translates to the 
upper bound of the constant term in the T = limit 
of Xkagome a.s c — 1.4 X 10"'^ cmu/mol-f.u. We fitted 
the SQUID data to Xbuik — c + Xcw below 3 K, and 
subtracted xhw ^o^i Xbuik to estimate Xkagome the 
entire temperature range. The obtained results in Fig. 2b 
are similar to the previous two cases, except for the finite 
intercept at T = 0. 

To summarize, we identified the source of a large en- 
hancement of spin susceptibility in ZnCu3(OH)6Cl2 as 
weakly interacting defect moments induced by Cu^^ in 
Zn sites and their immediate neighbors, and paved an 
avenue to analyze the bulk spin susceptibility data. By 
subtracting the defect contributions, we estimated an up- 
per and lower bound of the intrinsic spin susceptibility 
within the kagome planes, Xkagome (with a finite inter- 
cept) and Xkagome (with a small gap), respectively. The 
reality must be somewhere in between. 
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